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Huge Spin-Driven Polarizations at Room-Temperature in bulk BiFe0 3 


Jun Hee Lee* and Randy S. Fishman 

Materials Science and Technology Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, C/aSTQ 

The spin-driven polarizations of type-/ multiferroics are veiled by the preexisting ferroelectric 
(FE) polarization. Using first-principles calculations combined with a spin model, we uncover two 
hidden but huge spin-driven polarizations in the room-temperature multiferroic BiFeCU. One is 
associated with the global inversion symmetry broken by a FE distortion and the other is associated 
with the local inversion symmetry broken by an antiferrodistortive octahedra rotation. Compari¬ 
son with recent neutron scatterings reveals that first polarization reaches ~3.0 /iC/cm 2 , which is 
larger than in any other multiferroic material. Our exhaustive study paves a way to uncover the 
various magnetoelectric couplings that generate hidden spin-driven polarizations in other type-/ 
multiferroics. 

PACS numbers: 75.25.-j, 75.30.Ds, 75.50.Ee, 78.30.-j 


Although BiFe0 3 is endowed with a high ferroelec¬ 
tric (FE) and antiferromagnetic (AF) transition temper¬ 
atures, T c « 1100 K [I] and T N « 640 K |2], the disparity 
between T c and Tn T c in this type-/ multiferroic sug¬ 
gests that the magnetoelectric (ME) couplings may be 
quite weak. Despite enormous effort HE] , a microscopic 
picture embracing all of the ME coupling mechanisms 
in bulk BiFe0 3 is still missing. By contrast with type- 
II multiferroics where Tn = T c and the ME polariza¬ 
tions have been well characterized [7], the large FE po¬ 
larization, high Neel temperature, and long 62 nm period 
of BiFe0 3 have hindered measurement of its spin-driven 
ME polarization. Based on elastic m i and inelastic 
neutron-scattering [9], Raman-scattering [10], and THz 
spectroscopy m measurements of recently-available sin¬ 
gle crystals, it is now possible to provide detailed infor¬ 
mation about the intrinsic ME couplings in bulk BiFe0 3 . 
These results are crucial to control the electrical proper¬ 
ties of BiFe0 3 with a magnetic field and vice versa. 

Combining a first-principles approach with a spin- 
cycloid model, we explain the origin of all possible 
ME couplings and spin-driven (SD) polarizations pro¬ 
duced by exchange-striction (ES), spin-current (SC), and 
single-ion anisotropy (SIA). All polarizations are fostered 
by broken inversion symmetries with two types of lat¬ 
tice distortion in R3c bulk BiFe0 3 : FE and antifer¬ 
rodistortive (AFD). By comparing our results for the 
spin-driven atomic displacements with elastic neutron¬ 
scattering measurements 00112], we demonstrate that 
the ES-polarization (ESP) ^3 fiC/cm 2 dominates over 
other sources of polarization in the spin cycloid and is 
larger than any previously reported SD polarization. 

In type-/ multiferroics, the absence of an inversion cen¬ 
ter due to the preexisting FE polarization fosters the 


spin-driven polarizatioins. Specifically, the change of the 
scalar product • S j at the magnetic transition modu¬ 
lates the degree of broken-inversion symmetry and pro¬ 
duces corresponding ESPs m- While the FE distortion 
eliminates a global-inversion center, the AFD distortion 
eliminates a local-inversion center. Therefore, FE and 
AF distortions each generate their own ESP. 

All possible polarizations are obtained by differentiat¬ 
ing the Hamiltonian with respect to an electric field. For 
symmetric exchange couplings, ES is dominated by the 
response of the nearest-neighbor interaction Ji from the 
original Hamiltonian: 

n EX = -J2 S; • S ( = - ]T jf Si • s j, (1) 

(i,j) Ri ,Rj =Rj+efc 

where k = x, y, or z. Taking the FE polarization 
along z / =[111], the ESPs are then obtained from P ES = 
-V S U EX /N with 

P E ! =z'(z'-P ES ) 

= z'(2C± + C'n )z' • Wi = z'v^Cfe z' • Wi, (2) 

pES / p ES 

■* AFD — z x r 

= (C|—Ci)z'xW2= Cafd z'xW 2 , (3) 

Wik = 2 £ Si-Sj, (4) 

R-i )R-j —R-i -|-6fc 

W 2 k = L X (-l) ni Si-Sj, (5) 

R^ ,Rj = R.y 

where C± = 8Ji/8E a ((3 ^ a), and C\\ = 8Ji/8E a 
for spin bonds perpendicular and parallel to the electric 
field, respectively. 
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(a) G-AFM—»P = 93 |j,C/cm 2 (b) FM —•P = 99 j^C/cm 2 



FIG. 1: First-principles relaxation results for magnetic or¬ 
derings (a) C-AFM and (b) FM to calculate ESPs from FE 
distortion of Fe (z'(Fe)) and of Bi (z / (Bi))along [111] and 
AFD R 4 rotations ( 6 ) perpendicular to [111]. (0) and (3) are 
indices of hexagonal layers (x + y + z = n ). 


Unlike W\k , W%k alternates in sign due to opposite 
AFD rotations on adjacent hexagonal layers labeled by 
rii. The ESP parallel to z' with coefficient Cfe = 
(2C_l + C||)/a/3 modulates the FE polarization that al¬ 
ready breaks inversion symmetry above Tb; the ESP per¬ 
pendicular to z' has coefficient Cafd = Cj_ — Cy. The 
AFD breaks the local inversion between nearest-neighbor 
spins perpendicular to z' because each oxygen moves 
along [ 0 , 1 , 1 ], [ 1 , 0 , 1 ], and [ 1 , 1 , 0 ], perpendicular to z 
By ignoring the cycloidal harmonics but including the 
tilt [14| r produced by Pafd,^ the spin state propagat¬ 
ing along one of the three hexagonal orientaions x' can 
be approximated m as 

= S'( — l) n " +1 cost sin [2V2 tt 5 x.' -Hi/a ], (6) 
S y f(Iti) = S sin r sin [2v / 27r<5x'-R i /4 (7) 

S z >(Ri) = S'(-l)^ +1 cos [2>/27nfx^R*/a], (8) 

where a = 3.96 A is the pseudo-cubic lattice constant 
and S/(y/2a) = 62 nm is the period of the cycloid so that 
S ~ 0.0045. Recall that [Iff sinr = So/S where Mo = 
2/z b S 0 is the weak FM moment of the AF phase along y' 
above H c . For a moment [3] M 0 = 0.03//B, t = 0.006 or 
0.34°. Because higher harmonics are neglected, averages 
taken with the tilted cycloid introduce a very small error 
of order C 3 2 ~ 2.5 x 10 -5 . 

To calculate the electric-field induced lattice distortion 
and the associated change in J \, we employ a method PH 
that was successfully applied to other multiferroic oxides. 
The resulting ESP coefficients are given in Tab. [TJ Using 
the result Cfe = 215 nC/cm 2 and Eq.Q for Pf|, we find 
that the ESP for the simple twisted cycloid in Eqs.([ 6 ]- 


TABLE I: Calculated (LSDA+U) ESP components perpen¬ 
dicular C _l or parallel C\\ to the electric field and associated 
FE (Cfe) and AFD (Cafd) spin-driven polarization. Values 
in the parenthesis are directly obtained from the polarization 
difference of the relaxed structures with C-AFM or FM or¬ 
dering as shown in Eq. [Tol p~2|and Fig. [l] 


nC/cm 2 

C± 

c w 

Cfe 

Cafd 

LSDA+U 

186 

0.769 

215 (480) 

-185 (-108) 


[ 8 | is 

(Pf!) “ — Cfe5' 2 z / cos 2 r = —1.3 /iC/cm 2 z / . (9) 


Because a harmonic approximation m was used to gen¬ 
erate the possible polar distortions induced by the electric 
field, this ESP was only evaluated to quadratic order in 
the lattice distortions driven by the spin ordering. 

However, the ESP may be large enough to induce 
atomic displacements and lattice distortions beyond the 
harmonic limit. As shown in Fig. [TJ one can calculate 
the ESP more accurately including anharmonic effects 
and spin-lattice couplings by fully relaxing the atoms and 
the lattice for different magnetic orderings (G-AFM and 
FM) with 


A P T 


iES 

FE 


= Pfm — Pafm = 2Cfe£' 2 

= 7T ^2 — Z*AFM U i,AFM} 


i=Fe,Bi 


= 6.0 11 C/cm 2 , 


( 10 ) 


so that Cfe = 480 nC/cm 2 and 

(Ppi) = — Cfe*S' 2 z / cos 2 r = —3.0 /iC/cm 2 z / , ( 11 ) 


where Z*, and V represent Born effective charge, 
atomic position and volume, respectively. The change 
of spin ordering from G-AFM to FM shifts the Fe 
and Bi atoms by 0.020 A and 0.019 A, respectively as 
shown in Fig. [l] While the Bi effective charge hardly 
changes (Z^ FM =4.82e, Z FM =4.83e), the Fe effective 
charge changes significantly (Z^ FM =3.91e, Z FM =4.11e) 
due to the spin-induced hybridization between Fe and 
oxygen. Consequently, Cfe =480 nC/cm 2 is a factor of 
two larger than the harmonic value in Tab. |T| This result 
will later be compared with neutron-scattering measure¬ 
ments. 

The Supplement shows that Eq. © describes the ESP 
perpendicular to the FE-polarization direction z' due to 
the inversion-symmetry breaking from AFD rotations. 
Hence, the second ESP reflects the change of AFD and 
the associated local polarization driven by spins per¬ 
pendicular to zFor a simple tilted cycloid, (fU 2/e ) = 
(P/vfd) = 0 because the AFD distortions do not globally 
break inversion symmetry. 
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FIG. 2: Separation of FE (Dfe,/c) and AFD (Dafda) DM interactions in BiFeOs- Spin directions (red arrows) and DM vectors 
(blue arrows) are indicated, (a) and (b) have clockwise and counter-clockwise rotation of spins and use an 80 atom unit-cell 
(four spins in the unit-cell along each direction), (c) and (d) have zigzag-type spin ordering. +/— above red arrows denote 
signs of DM interaction and of spin-current between two nearest neighbor spins. 


As shown in Fig. [l] first-principles calculations can 
capture P E | D by evaluating the AFD-induced oxygen 
displacements perpendicular to z' with the change of 
spin ordering. The increase of the AFD rotation (A 0 = 
0.54°) from G-AFM (0 = 13.08°) to FM (0 = 13.62°) 
corresponds to an increase of the oxygen displacement 
(0.015 A) along [ 0 , — 1 , 1 ], [ 1 , 0 , — 1 ], and [— 1 , 1 , 0 ] perpen¬ 
dicular to zTherefore, 

APafd = Pfm - Pafm = 2 CafdS 2 

= ~ 3 ^° e 0.015l = -1340 nC/cm 2 , ( 12 ) 

so that Cafd = —108 nC/cm 2 . 

Interestingly, the two ESPs Cfe and Cafd are cou¬ 
pled. While the ESP components C± and C\\ coopera¬ 
tively increase ESP along z' under the inversion symme¬ 
try broken by FE, they produce opposite contributions to 
the AFD-induced ESP perpendicular to z'. Cj_ is largely 
positive due to the reduction of the Fe-O-Fe bond angle 
driven by the FE distortion, which favors FM from the 
Goodenough-Kanamori (GK) rules [18]; C\\ is almost zero 
because the bond contraction between Fe-O-Fe does not 
significantly alter the spin-density environment around 
the d 5 electrons of Fe. The large difference between C± 
and C || induces the large ESP induced by AFD rotations. 

The global and local inversion symmetry breaking by 
FE and AFD distortions produce the DM interactions 
DpE,k and Tafd,/c- Due to their distinct translational 
characters, they can be separated using the procedure 
sketched in Fig. [2] Since the FE distortion is globally 
uniform, its associated Dfe,/c is uniform too. Because the 
AFD rotation alternates between hexagonal layers, the 
associated DM vector Dafd,/c also alternates as shown 
by the blue arrows in Fig. [ 2 ] As shown in Figs. [ 2 ja) and 
(b), a cycloid consisting of four spins along e/~ generates a 
translation-invariant spin current x S j. The uniform 
component Tfe,/c is extracted from 

-®fe ]k = ^oi 4Dp E /c 5' 2 — - KS 2 , (13) 


^FE,k ~ gg2 (^FE ),k ^FE,k )> (^) 

where =L refer to counterclockwise (+) and clockwise (—) 
rotations, respectively. The translation-odd T)afd,/c does 
not appear in this expression. 

Using a zig-zag type spin arrangement that gener¬ 
ates a spin current S* x S j with alternating sign, the 
translation-odd -Dafd,/c is extracted from 

^AFD,fc = ^oi dD a FDj £ S 2 — -KS 2 , (15) 

^AFD,fc = g^2 (^AFD ,k ~ ^AFD ,fc)> 0-®) 

which does not contain the translation-even -Dfe,a- 
As for the ESP, the SC polarization (SCP) P sc = 
Pp E + Pafd splits into terms produced by the inversion- 
symmetry breaking of FE and AFD distortions, respec¬ 
tively. The first SCP is caused by the response of the FE 
distortion to an electric field: 


oo/SC 
tdSC _ un FE 

FE ’ 7 dE 1 

1 ®Dfe, k 

N ^ dE~ 

&k 1 


(Sj x Sj_|_ efe ). 


(17) 


The P-field derivatives of the DM, f k @ = dDvE,k/dEp 
are presented in Tab. [IT] and the Supplement 

The second SCP arises from the AFD rotation. Its sign 
alternates due to the alternating AFD rotation direction 
along [111]: 


oo/sc 

tdSC _ ori AFF> 

AFD ’ 7 dE 1 

( 18 ) 

e fc 7 


The SCP components a. k @ = dD afd, k/dEp are evalu¬ 
ated in the Supplement and presented in Tab. [II] 
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FIG. 3: (a) Spin-driven polarizations along z and AFD 

rotations perpendicular to z'. (b) Calculated Fe polarizations 
along z induced by ES (P ES ) compared with measurements 
by Lee et al. [6], Arnold et al. [5], and Palewicz et al. H2. 
(c) Bi polarizations along z! induced by ES. (d) Total SD 
polarization A P z > (F e)+AP z / (Bi). Inset is the SC polarization 
calculated from Eqs. ( |21|22| and Ref. [20] . (e) Change of AFD- 
rotation angle induced by ES. (f) AFM spin ordering versus 
temperature calculated by mean-field theory. 


SIA is the last possible source for the ME polarization. 
Starting with the SIA enenrgy P SIA = — K JA(S^ • z') 2 , 
the SD-polarization has z' component 

/ pSiA \/3£ „'\2 c hm 

zP =-r 2 A s *'*> • (=■$£:■ (19) 

i 


TABLE II: Calculated (LSDA+P) SCPs and SIA polariza¬ 
tion (SIAP) derived from DM and K interactions. Subscripts 
denote whether the spin-bond direction is parallel or perpen¬ 
dicular to the electric field. f xy , f xy are not shown due to 
low convergence. 





SCP from D 



SIAP 

nC/cm 2 


«r 

a y V 

a x x y 

a x z y 

nxy 

Jy 

4 

LSDA+P 

6.6 

-21 

16 

1.2 

-13 

-4.9 

16 


The elongation £ of K is given in Tab. [Hi 

For the simple tilted cycloid of Eqs.(|6][8]), 

< pSIA> = 7 §ZS 2 z\ (20) 

(Pfe) = —27 tV6S 2 5cost 

{(/r-/D y' + ^(fz V ~fy V ~fy X ) Z '}> (2!) 

(n xx 1 

(P| ¥D )=-^ S i^{ i - + „ r+ a ?+ „.. + <«-},'.( 22 ) 

Symmetry relations for and a%P are given in the 
Supplement. Because S, r <C 1, the projected polariza¬ 
tion along z' is z' • (P sc + P SIA ) w z' • (P SIA ) w 87 
nC/cm 2 which is larger than the experimental value 40 
nC/cm 2 [3] 0] obtained from the jump in polarization 
below the critical field H c ^ 20 T. Intriguingly, RpE.fc 
produces an additional SCP along y', which may explain 
why the SC also generates a polarization perpendicular 
to z' [4]. Since the SCP and SIAP are still much smaller 
than the ESP, the dominant polarization at the magnetic 
transition is driven by ES. 

Figure [3] shows all the ME polarizations driven by 
the AFM spin-ordering around Xn and compares those 
results to elastic neutron-scattering measurements [6j 
0 [12]. Although the neutron-scattering data is rather 
spread, all three papers indicate that both the polariza¬ 
tion and AFD rotation angle are reduced by the huge ES 
around Tn- As explained in the Supplement, we convert 
the preliminary neutron-scattering data to the change 
of the spin-driven polarization at Xn using Ginzburg- 
Landau free energies. With spin ordering (Fig. (3^d)), Fe 
and Bi move —0.010A and —0.0095A, respectively, and 
induce polarizations AP(Fe)= —1.7 /iC/cm 2 and AP(Bi) 
= —1.3 /iC/cm 2 . The net induced polarizaton P tot = 
AP(Fe)+AP(Bi) = —3.0 /iC/cm 2 is in excellent agree¬ 
ment with neutron-scattering measurements. 

Intriguingly, both Bi and Fe shift with the spin order¬ 
ing. Fe moves antiparallel to reduce the FE polariza¬ 
tions because AFM ordering favors 180° displacements 
by the GK rules m ■ The reduction in the Fe po¬ 
larizations simultaneously reduces the Bi polarization. 
Consequently, the magnitude of the net polarization is 
greater than any previously-reported spin-driven polar¬ 
ization (0.29 /iC/cm 2 in CaMn 70 i 2 [19] and 0.36 /iC/cm 2 
in GdMn 205 ED])- 
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We have also discovered another huge but hidden 
ES due to AFD rotations that are strongly coupled 
to the FE-driven ES. Obviously, this contribution can¬ 
not be easily measured because AFD rotations do not 
break global inversion symmetry and do not produce a 
net macroscopic polarization. However, AFD rotations 
do break local inversion symmetry and their associated 
atomic displacements appear in the neutron-scattering 
data in Fig[3^e). 

Based on the good agreement between our predic¬ 
tions and neutron-scattering results, we conclude that the 
AFD rotation angle is suppressed by ES. Both the polar¬ 
ization and AFD reduction can be understood in terms of 
the GK rules [18]: AFM ordering decreases bond angles, 
which reduces the FE polarizaton and AFD rotations. 
Due to recent advancement in local polarization mea¬ 
surements m , it may soon be possible to directly image 
the spin-driven structural modification of the AFD. 

Although some calculations predict polarizations ^ 
6 /iC/cm 2 {22, 23 for orthorhombic perovskites such as 
HoMnOs, the largest measured spin-driven polarization 
prior to this work was found in pressurized TbMnOs, 
where P = 1.0 /iC/cm 2 can rise to 1.8 /iC/cm 2 with 5.2 
GPa [24] at 5 K. We have checked the possible spin-driven 
polarization in another type-/ multiferroic, BiCoOs, from 
neutron scattering [25] . But its polarization appears to 
be smaller than that of BiFeOs by one order of magni¬ 
tude. Hence, the type-/ multiferroic BiFeOs unexpect¬ 
edly exhibits the largest ever spin-driven polarization 
3 /i C/cm 2 ) at room temperature. There are three 
reasons for this huge spin-driven polarization. First, the 
exchange interactions and their response to external per¬ 
turbations such as electric field or temperature are larger 
than for the DM interactions, as shown in Tab. [T] and 
Tab.|H| Second, even if the ES coefficients were the same 
for BiFeOs and the above manganites, the enhanced spin 
correlation function (|(Si • Sj) |=P 2 =6.25 in BiFeOs , and 
2.25 in the manganites) would strongly enhance the ESP 
in BiFeOs. Third, in contrast to the P-type ordering 
in the manganites, the almost antiparallel alignment of 
neighboring spins in BiFeOs also enhances the ESP. 

The greatest advantages of BiFeOs are its large FE po¬ 
larization, high T c , and high Tn above room temperature. 
These advantages have unfortunately hampered precise 
characterizations of the ME polarizations around Tn- 
Leakage currents at high temperatures and the preexist¬ 
ing large FE polarization have hidden the spin-driven ME 
polarizations at Tn- Fortunately, intrinsic measurments 
such as neutron-scattering, Raman spectroscopy and di¬ 
rectional dichroism have recently begun uncovering the 
hidden ME couplings of BiFeOs. So in addition to having 
the largest known FE polarization, BiFeOs ma y also have 
spin-driven polarizations much larger than in any other 
known material. Our systematic approach will greatly 
aid further exploration of hidden but possibly large spin- 
driven polarizations and their ME origins in other type-/ 


multiferroics. 
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